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Abstract  projectile's  aerodynamic  stability  will  be  consi¬ 

dered.  Within  that  context,  the  flare  angle  and 
The  effect  of  nose  shape  on  two  flare  the  flare  length  will  be  varied  to  study  the  effect 

stabilized  projectiles  was  studied  using  a  of  these  two  parameters.  A  brief  description  of 

Parabolized  Na  v  i  e  r  -  S  t  ok  e  s  code.  Pressure  the  PNS  code,  code  inputs,  configurations,  and 

coefficients,  forces  and  moments ,  skin  friction  flow  conditions  will  be  presented  in  the  second 

coefficients,  and  Stanton  number  calculations  are  section.  Aerodynamic  results  will  be  compared  for 

presented  for  the  hemi sphere -cy  1  i nder- f  1  are  and  each  configuration  in  the  third  section.  Finally, 

the  cone-cylinder-flare  configurations.  Pitching  conclusions  to  be  drawn  will  be  presented  in  the 

moment  and  static  margin  plots  versus  flare  angles  last  section, 

and  flare  lengths  are  presented  in  a  parametric 

study  to  show  aerodynamic  stability  effects.  PNS  Code  and  Configurations 

Distinctly  different  flow  field  values  of 

pressure,  local  Mach  number,  and  dynamic  pressure  A  Parabolized  Na v ier-S t okes  code”  used  for 

were  generated  by  the  two  different  nose  this  numerical  study  generates  finite  difference 

configurations.  These  flow  field  values  just  approximations  to  steady,  three-dimensional 

upstream  of  the  flare  will  be  examined.  The  code  solutions  of  supersonic  flow  over  arbitrarily 

demonstrated  its  value  as  a  design  tool  by  making  shaped  bodies  at  high  Reynolds  number  as  long  as 

a  clear  distinction  between  aerodynamics  1 1  y  stable  there  are  no  large  subsonic  or  axially  separated 

characteristics  for  this  variety  of  nose  shapes,  regions.  From  an  initial  condition  starting  solu- 


flare  angles,  and  flare  lengths. 

Introduction 

The  stabilizing  effect  of  flared  afterbodies 
has  long  been  incorporated  in  the  design  of 
rockets  and  projectiles.  For  short  projectiles  at 
relatively  high  Mach  numbers,  the  flare  angle  and 
the  flaie  length  can  play  a  significant  role  in 
the  projectile's  aerodynamic  stability.  Extensive 
experimental  studies  have  been  made  on  cone- 
cylinder-flared  projectiles  at  Mach  numbers 
ranging  from  2  to  4  1  and  a  Parabolized  Navier- 
Stokes  (PNS)  code  has  been  used  for  numerical 
studies.^*'’  Most  of  these  studies  involved 
conical  nosed  bodies  at  lower  supersonic  speeds. 

This  paper  will  present  aerodynamic 
characteristics  determined  numerically  for  two 
flare  stabilized  conf igurations  at  Mach  7  and  at 
2”  angle  of  attack;  a  hemisphere-cylinder-flare 
and  a  cone -c y 1 i nder - f 1  are.  Mach  7  was  chosen 
because  it  was  high  enough  to  yield  characteris- 
ticly  high  Mach  number  effects  yet  low  enough 
that  real  gas  effects  can  be  ignored.  The  low 
angle  of  attack  was  selected  to  generate  pitching 
moment  data.  The  effect  of  the  nose  shape  on  each 
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tion,  the  code  marches  downstream  from  any  given 
axial  location  on  the  body.  For  the  blunt  hemis¬ 


pherical  nose  configuration  in  Figure  1(a),  the 
code  needed  an  input  starting  solution  obtained 
from  a  thin-layer  unsteady  Navier-Stokes  code. 


For  the  sharp  conical  nose  configuration  in  Figure 
1(b),  the  PNS  code  has  a  self-start  capability 
wherein  initial  starting  planes  are  calculated 
within  the  code. 


The  PNS  code  is  a  space  marching  code  which 
uses  a  numerical  method  originally  developed  at 
NASA  Ames  Research  Center  by  Schiff  and  Steger. 
This  technique  obtains  finite  difference  solutions 
by  solving  implicit  difference  equations 
noni tera t i vely  by  way  of  local  linearization  of 
the  flux  vector.  The  resulting  system  of 
algebraic  equations  are  solved  using  the  Beam- 
Warming  method^  of  implementing  the  alternating 
direction  implicit  scheme  using  the  delta 
formulation.  Several  additional  modifications  to 
the  code  have  been  made.  J 


The  PNS  flow  field  predictions  for  the 
hemi s phe re - c y 1 i n d e r - f  1  a  r  e  were  obtained  by 
marching  from  a  starting  solution  generated  by  a 
thin-layer  unsteady  Navier-Stokes  code.  This 
starting  solution  was  obtained  at  an  axial 
location  at  the  hemisphere-cylinder  shoulder  with 
an  algebraic  grid  of  19  by  30  points  in  the 
circumferential  and  radial  directions, 
respectively.  This  data  was  then  interpolated  to 
obtain  an  algebraic  grid  of  19  by  45  points.  The 
data  was  obtained  at  0°  angle  of  attack  and  then 
rotated  to  a  2°  angle  of  attack.  For  a 
hemispherical  nose,  this  rotation  presents  no 
problems  theoretically.  Thirty-five  of  the  45 
radial  grid  points  were  clustered  near  the  wall  to 
resolve  the  boundary  layer.  A  marching  step  size 
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of  0.01  cylinder  diameters  W3S  used  up  ro  a  loca¬ 
tion  one-half  cylinder  diameter  aft  of  the 
hemi s phe re -cy 1 i nde r  shoulder.  This  small  step 
size  captured  the  strong  gradients  in  the  flow 
field  in  this  shoulder  region.  The  step  size  was 
then  increased  to  0.03  diameters  for  marching 
along  the  cylinder  to  a  plane  0. !  3  diameters 
upstream  of  the  cy  1  i nde r- f 1  are  juncture.  A 
reduced  step  size  of  0.02  diameters  was  used  to 
adjust  for  the  rapid  pressure  rise  in  the  region 
at  the  flare.  This  step  size  was  close  to  one- 
third  the  boundary  layer  thickness  predicted  by 
the  PNS  cod t*  at  this  point.  Analysis  by  Vi.  gas 
and  Horst  man  indicate  that  a  step  size  one-half 
the  boundary  layer  thickness  is  sufficient  to 
resolve  two-dimensional  shock /bou ndary  layer 
interactions.  Calculations  were  made  on  0,  5, 
10,  15,  20,  and  2  5  degree  flare  angles.  Attempts 
to  calculate  a  30"  flare  angle  failed  with  axial 
flow  separation  characteristics.  This  is 
consistent  with  the  PNS  code's  ability  to  detect 
the  onset  of  incipient  separation  as  evaluated  by 
Cottrel  1  and  Chapman.^  This  failure  may  also  be 
attributed  to  the  Beam-Warming  type  of  algorithm 
which  tends  to  introduce  errors  in  the  form  of 
local  flow  property  oscillations  as  a  result  of 
the  cent ra l -d i f f erenc ing  of  fluxes  across  discon¬ 
tinuities.^  Versions  of  the  PNS  code  using  an 
upwind -differencing  scheme  to  alleviate  this 
problem  are  being  employed.  The  code's  smoothing 
terms  were  as  follows:  explicit  damping  parameter 
of  0.2,  implicit  damping  parameter  of  0.4,  sta¬ 
bility  coefficients  of  0.2  5.  The  nominal  value  of 
the  grid  spacing  at  the  wall  was  0.5  x  10”  . 

The  PNS  flow  field  prediction  for  the  cone- 
cy  l  i nder- f l are  configuration  was  obtained  by 
marching  from  a  solution  generated  by  the  code's 
self-start  option  for  sharp  conical  noses.  This 
starting  solution  was  obtained  at  an  axial 
location  1  cylinder  diameter  from  the  nosetip  at  a 
2  °  angle  of  attack.  The  algebraic  g r i d  w a s 
19  by  45  points  in  the  circumferential  and  radial 
directions,  respectively.  Thirty-five  of  the  45 
radial  grid  points  were  clustered  near  the  wall  to 
resolve  the  boundary  layer.  This  grid  is  set  up 
exactly  the  same  as  the  grid  for  the  hemi sphe re  - 
cylinder-flare.  However,  the  volume  of  the  grid 
for  the  cone -cy l  i  nder- f 1  are  case  is  much  smaller 
than  the  volume  of  the  grid  for  the  hemi sphe re- 
cylinder-flare  case  due  to  the  smaller 
distance  of  the  bow  shock.  Hence,  the  grid  points 
clustered  near  the  wall  for  the  cone-cv  l  inder- 
flare  are  packed  closer  together.  A  marching  step 
size  of  0.02  diameters  was  used  to  a  plane  0.8 
cylinder  diameters  aft  of  the  cone-cylinder 
shoulder.  This  small  step  size  resolved  the  large 
gradients  in  the  flow  field  as  the  solution 
marched  up  over  the  shoulder.  The  step  size  was 
then  increased  to  0.04  diameters  for  marching 
along  the  cylinder  to  a  plane  0.14  diameters 
upstream  of  the  cylinder-flare  juncture.  The  step 
size  was  reduced  back  down  to  0.02  diameters  to 
adequately  resolve  the  rapid  pressure  rise  in  the 
flare  region.  For  the  20*  flare,  the  step  size 
was  kept  at  0.04  diameters.  A  flare  angle  of  25° 
was  attempted  but  the  code  indicated  the  onset  of 
axial  flow  separation  just  as  the  30*  flare  angle 
attempt  for  the  hemisphere-cylinder-flare  did  and 
the  reasons  cited  at  that  time  are  applicable 
here,  also.  The  smoothing  terms  were  kept  the 
same  as  those  used  in  the  hemi sphe re -cy  1  i nder - 
flare  study  for  all  cases. 


Results  and  Discussion 


He mi sphe re -Cy linder-Flare 

Figure  2  shows  the  shock  structure  for  the 
hemisphere-cylinder-15  degree  flare  at  Mach  7. 
The  outer  boundary  shock  was  fit  by  the  unsteady 
Na vier-S tokes  code  for  the  hemisphere  nose  segment 
and  the  PNS  code  fit  the  shock  for  the  cylinder- 
flare  segment.  Note  the  small  discontinuity  where 
the  shock  fits  are  joined  together  at  the  shoulder 
region.  This  is  not  unusual  when  shifting  f rom  an 
unsteady  Na  v  i  e  r -S  t  ok  es  code  to  a  PNS  code.  The 
unsteady  code  captured  the  subsonic  region  of  the 
nose.  The  sonic  lines  were  estimated  from 
locations  where  the  total  Mach  equaled  1.  The 
embedded  secondary  shock  locations  were  captured 
by  the  PNS  code  and  the  estimated  locations  in  the 
flare  region  were  obtained  from  P / P •  f 
distributions  between  the  body  and  the  shock  at 
the  leeward  and  at  the  windward  planes.  These 
captured  shocks  are  smeared  across  grid  points  and 
sharp  discontinuities  are  not  apparent.  Figure  3 
shows  pressure  coefficient  distributions  at  the 
surface  at  three  circumferential  stations:  a 
windward  ray,  a  90°  or  side  ray,  and  a  leeward 
ray.  As  anticipated  at  this  sma  l  L  angle  of 
attack,  there  is  a  slightly  larger  Cp  on  the 
windward  side  which  decreases  slowly  as  one  moves 
circumferentially  to  the  lower  Cp  on  the  leeward 
side.  From  a  Cp  of  1.83  at  the  nose  stagnation 
point,  the  flow  field  pressure  undergoes  a  rapid 


circumferentially  to  the  lower  Cp  on  the  leeward 
side.  From  a  Cp  of  1.83  at  the  nose  stagnation 
point,  the  flow  field  pressure  undergoes  a  rapid 
expansion  over  the  nose  shoulder,  decreases 
gradually  along  the  cylinder  body,  and  jumps 
quickly  at  the  cy  1  i  nde  r -f  1  a  re  juncture.  Not**  once 
again  that  the  joining  of  unsteady  Na v ier-Stokes 
code  data  with  PNS  code  data  results  in  a  slight 
Cp  discontinuity  in  the  shoulder  region.  Figures 
4-6  show  pressure  coefficient  distributions  at  the 
same  three  circumferential  stations  on  the 
cylinder-flare  segment  for  flare  angles  at  0  and 
15,  5  and  20,  and  10  and  2  5  degrees.  Note  how  the 
pressure  jumps  at  the  c v 1  i n d e r - f 1  a r e  juncture. 
For  lower  flare  angles,  the  pressure  jumps  and 
then  decays  as  expected  for  uniform  flow 
approaching  the  flare  (i.e.,  jump  to  pressure 
associated  with  a  2-D  wedge  and  then  drop  down  to 
pressure  associated  with  a  cone).  At  the  20  and 
25  degree  flare  angles,  there  appears  a  deviation 
from  this  trend.  This  is  due  to  nonuni  form  flow 
a  p p  r oa c h i n g  t  h  e  flare  and  will  he  discussed 
further  in  the  comparison  with  the  cone-cylinder 
case.  Attempts  to  calculate  a  30"  flare  cast' 
failed  with  axial  flow  separation  characteristics 
for  reasons  discussed  earlier. 

It  should  be  noted  that  this  marching  code 
cannot  predict  the  pressures  on  the  ease  of  tin* 
flare.  If  the  pressure  distribution  on  the  base 
is  uniform,  then  the  base  pressure  effect  on 
pitching  moment  is  negligible  for  this  small  angle 
of  attack.  Pitching  moment  curves  versus  flare 
lengths  for  various  flare  angles  are  shown  in 
Figure  7.  Stability  grows  with  increases  in  both 
flare  length  and  flare  angle.  These  results  are 
similar  to  the  stability  trends  reported  by 
Cottrell  and  Chapman. 

Static  margin  curves  were  constructed  by 
subtracting  the  center  of  gravity  (eg)  location 
from  t  h  e  center  of  pressure  location.  A 
homogeneous  b  o  d  y  was  assumed  for  the  eg 
calculation.  The  center  of  pressure  locations 
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were  calculated  by  the  PNS  code.  Both  locations 
were  measured  positively  from  the  nose.  When  the 
center  of  pressure  is  farther  away  from  the  nose 
than  the  eg,  the  nose  will  pitch  downward  which 
indicates  aerodynamic  stability.  In  Figure  8,  the 
static  margin  for  0,  5,  10,  15,  20,  and  25  degree 
flare  angles  is  plotted  for  various  flare  lengths. 
For  the  0,  5,  and  10  degree  flare  cases,  these 
curves  indicate  unstable  aerodynamics  for  whatever 
flare  length  chosen.  The  10°  flare  indicates  an 
impasse  where  an  increase  in  flare  length  moves 
the  eg  and  the  center  of  pressure  the  same 
distance  with  no  further  stability  deterioration 
or  increase.  Note  that  the  15,  20,  and  25  degree 
flare  cases  go  stable  rather  quickly  with 
increasing  flare  angle  and  flare  length  in  a 
manner  similar  to  that  reported  by  Cottrell  and 
Chapman .  Th  ‘r«£  three  flare  angles  all  go 
stable  for  fla*e  lengths  over  one-third  cylinder 
diameter  tong. 

Figures  9  and  10  show  typical  skin  friction 
coefficient  and  heat  transfer  (Stanton  number) 
distributions.  In  Figure  9,  skin  friction  decays 
smoothly  with  axial  position  except  near  the 
cy l  inder -f 1  are  juncture.  The  expansion  of  the 
flow  at  the  cone-cylinder  shoulder  increases  the 
boundary  layer  thickness  resulting  in  a  rapid 
decrease  in  skin  friction  coefficient  which 
steadies  to  a  smooth  decay  down  the  cylinder 
segment.  The  spike  in  the  skin  friction 
coefficient  at  the  cy 1 i nder - f  1  a r e  juncture  is 
physically  unreal.  It  may  be  attributed  to  the 
nature  of  the  PNS  type  of  algorithm. 
In  Figure  10,  Stanton  number,  or  heat  transfer, 
also  decreases  smoothly  with  axial  position  except 
near  the  cy  1  i nde r - f 1  a  re  juncture.  The  spike  at 
that  juncture  is  physically  unreal. 

Cone -Cylinder -FI are 

Figure  11  shows  the  shock  structure  for  the 
cone -cy l i nde r- 1  5  degree  flare  at  Mach  7.  The 
attached  shock  at  the  nose  and  the  outer  shock 
location  were  estimated  by  the  PNS  code  in  the 
step-back  mode  up  to  a  location  one  cylinder 
diameter  from  the  nosetip.  Planes  of  data  were 
stored  there.  The  PNS  code  picked  this  data  up 
and  marched  down  the  rest  of  the  configuration. 
Once  again,  the  estimated  embedded  shock  locations 
on  the  flare  were  approximated  from  the  P / P  £  n  f 
distributions  between  the  body  and  shock  at  the 
windward  and  leeward  planes.  Figure  12  shows 
pressure  coefficient  distributions  along  the 
configuration  at  the  surface  for  the  15°  flare 
case.  Note  the  Prandt  l  -  Meyer  expansion  at  the 
cone -cy 1 i nder  shoulder  which  drops  pressure 
rapidly  from  a  stagnation  Cp  of  2.13. 
Then  there  is  a  gradual  expansion  back  to  near 
freestream  pressure  values  down  the  cylinder 
segment.  Figures  13-15  show  pressure  coefficient 
distributions  on  the  cylinder-flare  segment  for 
flare  angles  of  0  and  15,  5  and  20,  and  10 
degrees.  Note  how  the  pressure  data  jumps  at  the 
cylinder-flare  juncture  for  all  cases  followed  by 
a  continual  pressure  rise  normally  associated  with 
nonuniform  flow.  This  nonuniform  aspect  of  the 
flow  will  be  discussed  in  the  comparison  section 
that  follows.  Attempts  to  calculate  a  25°  flare 
failed  with  axial  flow  separation  cha racter i s t  ics 
consistent  with  the  results  found  by  Cottrell  and 
Chapman 4  but  may  also  be  attributed  to  the  Beam- 
Warming  algorithm  mentioned  previously.  Pitching 


moment  curves  versus  flare  lengths  for  various 
flare  angles  are  shown  in  Figure  16.  As  expected, 
the  pitching  moment  grows  as  flare  angle  and  flare 
length  increases. 

Static  margin  curves  were  constructed  from  the 
PNS  code  data  and  are  shown  for  flare  angles  0,  5, 
10,  15,  and  20  degrees  in  Figure  17.  Note  it 
appears  that  the  0  and  5  degree  flare  angles  will 
never  yield  aerodynamic  stability  by  increasing 
the  flare  length.  Note,  too,  that  the  10°  flare 
goes  stable  at  a  flare  length  approxi ma t e 1 y  one- 
third  cylinder  diameter.  The  15  and  20  degree 
flares  yield  stability  rather  quickly  as  the  flare 
length  increases. 

Figures  18  and  19  show  skin  friction 
coefficient  and  heat  transfer  (Stanton  number) 
distributions.  The  spike  at  the  flare  for  both  of 
these  plots  is  not  realistic  and  may  be  attri¬ 
buted  to  flow  property  oscillation  problems 
mentioned  earlier. 

Comparison  of  Conf igurat ion  Resu Its 

A  comparison  of  the  aerodynamic  results  must 
start  with  an  examination  of  the  outer  shock 
boundaries  in  Figures  2  and  11.  The  hemisphere 
nose  generates  a  detached  bow  shock  with  a 
characteristic  subsonic  region  while  the  cone  nose 
has  an  attached  shock.  The  difference  between 
these  bow  shock  waves  lead  to  very  different  flow 
conditions  in  front  of  the  flares.  Figures  20, 
21,  and  22  show  the  radial  distribution  of  P/P£nfi 
local  Mach  number,  and  dynamic  pressure  (q  /q  ^ n ^ » 
respectively.  These  distribut ions  are  taken  at  a 
location  just  upstream  of  the  cylinder-flare 
juncture  for  the  hem i s phe re -cy 1 i nde r - f 1  a  re  and 
for  the  cone-cylinder-flare  configurations.  Also 
shown  in  these  figures  are  the  effective  height  of 
a  10°  and  a  20®  flare  at  the  base.  This  gives  a 
rough  idea  of  the  flow  encountered  by  these 
flares.  Figures  23(a)  and  23(b)  show  the 
algebraic  grid  at  these  locations  for  the 
hemisphere-cylinder-flare  and  for  the  c  one- 
cy  1  i  nder-f  1  are  ,  respectively. 

The  P/Pjnf  plot  in  Figure  20  reveals  the  large 
region  of  constant  pressure  each  configuration 
has.  Pressures  just  before  the  flare  are  close  to 
50%  higher  for  the  hemispherical  nose  configura¬ 
tion.  The  strong  bow  shocks  on  each  of  the  two 
lead  to  significant  entropy  losers.  The  blunt 
hemisphere  body  has  the  greater  losses.  This  is 
illustrated  by  the  lower  local  Mach  number  and  the 
lower  dynamic  pressure  in  front  of  the  flare  for 
this  blunt  nose  compared  to  the  cone  nose  as  shown 
in  Figures  21  and  22.  The  significant  reduction 
in  dynamic  pressure  accounts  for  the  much  lower 
static  margins  of  the  hem i s ph e r e -cy l inder -f  1  are 
compared  to  the  cone -cy l i nder- f 1  are  configuration. 

Although  the  hemispherical  nosed  body  has  the 
greater  entropy  losses  of  the  two,  the  conical 
nosed  body  has  more  of  a  variation  in  local  Mach 
number  radial  ly  over  the  distance  comparable  to 
the  flare  height.  It  is  this  variation  that 
accounts  for  the  rather  smooth  pressure  coeffi¬ 
cient  rise  along  the  flare  on  the  conical-nosed 
body  seen  in  Figures  12-15.  This  variation  in 
local  Mach  number  is  only  present  in  a  small 
radial  region  for  the  hemispherical -nose  body. 
Hence,  only  the  smaller  flare  angles  for  the 
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hemispher ica  1  -nose  hotly  have  smooth  prossuro  rises 
on  the  flare  as  seen  in  Figures  3-6. 

Cone lud ing  Remarks 

Based  on  the  predictions  made  with  the  PNS 
code  on  the  hemisphere-cy l i nder-f l are  at  various 
flare  angles  and  flare  lengths  with  flow 
conditions  of  Mach  7  and  at  2°  angle  of  attack, 
the  following  conclusions  can  bo  drawn. 

At  this  hypersonic  speed,  the  nose  shape  has  a 
distinct  effect  on  the  aerodynamic  stabilizing 
effect  of  flared  afterbodies.  Of  the  two  configu¬ 
rations  examined,  the  con i ca 1 -nosed  body  has 
greater  stability  and,  therefore,  more  utility  at 
hypersonic  applications.  It  appears  an  attached 
shock  wil  1  keep  entropy  losses  to  a  minimum.  This 
provides  the  stabilizing  control  surfaces  within 
tlie  outer  shock  enveloping  the  body  a  more 
effective  role  in  the  projectile's  hypersonic 
mission.  Furthermore,  the  increased  stability 
obtained  at  larger  flare  angles  and  greater  flare 
lengths  for  both  configurations  have  been 
p  resen  ted. 
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x  =  0.5d  /cn=12d 

4>f  =  0°  to  25°  /cy  =  2.13d 

*  =  22.50°  /  =  0.19,0.51,0.83d 

1 1 

Figure  1 .  (a)  Hemisphere  -  Cylinder  -  Flare  Geometry 
(b)  Cone  -  Cylinder  -  Flare  Geometry 
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Figure  2.  Shock  shape  of  hemiaphere-cy l inder- 
flare  for  15*  flare. 
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Figure  5.  Pressure  coefficient  distribution 
for  henisphere-cy lioder-f iare  for  5*  and  20* 
flares. 
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Figure  3.  Pressure  coefficient  distribution 
for  he»isphere-cy linder-f tare  for  15  flare. 
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Figure  6.  Pressure  coefficient  distribution 
for  he«isphere~cy 1 inder-f lare  for  10  and  25 
flares. 
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Figure  4.  Pressure  coefficient  distribution 
for  heaisphere-cy 1 inder-f Isre  for  0*  and  15* 
f I  ares . 
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Figure  7.  Pitching  sownt  curvea  for 
he.isphere-cy 1 inder-f 1  are  for  0,  5,  10,  15, 
20,  and  25  degree  flares. 
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Figure  8-  Static  margin  curves  for 
he.isphere-cylinder-flare  for  0,  5,  10,  15, 
20,  and  25  degree  flares. 
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Figure  9-  Coefficient  of  friction  curves  for 
hemisphere-cy linder-f lare  for  15*  flare. 
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Figure  10.  Stanton  number  (heat  transfer)  for 
hemisphere-cy linder-f lare  for  15*  flare. 
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Figure  11.  Shock  shape  for  cone-cy linder- 
flare  for  15*  flare. 
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Figure  12.  Pressure  coefficient  distribution 
for  cone-cylinder-flare  for  15*  flare. 
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Figure  13.  Pressure  coefficient  distribution 
for  cone-cylinder-flare  for  0*  and  15*  flares. 
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Figure  23.  Algebraic  grids  located  at  data 
planes. 
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at  the  windward  plane. 
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